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A ferromagnet-normal-metal-ferromagnet trilayer does not conserve spin current and as a 
3 , result, the conduction electrons can create a torque on the magnetic moments. When the 

trilayer is connected to a superconducting electrode, the spin torque created by a current 
Cn I flowing perpendicularly to the trilayer can drive the system to a configuration where the two 

magnetic moments are perpendicular to each other. Here we argue that in contrast to the non- 
equilibrium torque, the equilibrium torque (or magnetic exchange interaction) does not single 
r^ , out the perpendicular configuration. We calculate the equilibrium torque for a one-dimensional 

Q ■ lattice model and find that it is sensitive to the presence of the superconductivity. 
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1 Introduction 

Spin current is a useful concept to describe the effect of the conduction electrons on the dy- 
namics of the magnetic moments in a ferromagnet-normal-metal-ferromagnet trilayer. EJoBq The 
interaction with the local exchange field causes the potential felt by the conduction electrons to 
^ ' be spin dependent. In the presence of a spin-dependent potential, the spin current carried by 

Q ■ the conduction electrons is not conserved in the ferromagnetic layers. Since the total spin of the 

^ O^ i system is conserved, the missing spin current must be accounted for as a torque acting on the 

magnetic moments of the ferromagnetsp The torque from non-conservation of the equilibrium 
spin currents is known as the magnetic exchange (or RKKY) interactionp (There is a direct 
analogy between the equilibrium spin current flowing in presence of an non zero angle between 
\^ [ the two magnetic moments and the persistent current flowing through a mesoscopic ring in 

^^ ■ presence of an Aharonov-Bohm flux.) An additional torque is exerted out of equilibrium when 

(^ . a current is passed through the system. Experimentally, this non-equilibrium torque has been 

shown to be large enough to switch the relative orientation of the two magnetic moments of the 
two magnetic layers from parallel to antiparallel to each other. El In this paper, we report on 
jrt ■ our investigation of these torques when the trilayer is connected to a superconducting electrode. 

The non-equilibrium torque has been considered in Ref. |6| where it is shown that the torque 
can drive the system to a configuration where the two magnetic moments are perpendicular 
to each other. Here, we address the question whether a similar effect can be expected for the 

O . (equilibrium) magnetic exchange interaction. 
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k> ! 2 Definition of the spin torque 



Following Ref. ^, we consider a trilayer system consisting of two ferromagnetic layers Fq and F^ 
with a normal metal (N) spacer and one superconducting contact (S), see Fig. |I]. The magnetic 
moment of the layer F;, adjacent to S is considered fixed, e.g., by anisotropy forces. In a lattice 
formulation, the trilayer is described by the the Bogoliubov-De Gennes (BdG) equation^ 

Hi^i + tJ2'^j = e^i- (1) 



Here ^i is the quasiparticle wavefunction amplitude on site i = {ix,iy,iz), which has an electron- 



hole as well as a spin index, ^j = (■(/'i , ^j , V'i \i^i )'^- Further, t is the hopping amplitude, 
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Figure 1: Schematic of the Ferromagnet-Normal-metal-Ferromagnet (FNF) trilayer under consideration. Two 
ferromagnetic layers, Fa and F;,, separated by a normal metal spacer (N), are connected to one normal- metal 
contact (left), and one superconducting contact (S, right). (A normal-metal spacer is added between Fj and S.) 
The arrows indicate the directions of the local exchange fields, ha (6i = 6 for i £ Fa) and ht, (Oi = for i G F;,). 
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where Vi is an on-site potential, ep the Fermi energy, hi the local exchange field which we take in 
the xz plane, and 9i the angle between /i, and the z axis. The superconducting gap Aj = Aq in 
S and Aj = elsewhere. (This is appropriate if S is coupled to the trilayer by a point contact.) 
The spin current Jj^ is the current associated with the time derivative of the spin density. 
The spin current Jj^ g carried by a quasi-particle at energy e and with wavefunction ^ reads 



4,e = -tlmJ2 {^f+a^ri - ^'^Ml, 



(3) 



with d = (1, 0, 0), while the equilibrium spin current Jl^ is 
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e>0 



Here, a = {crx,<Ty,(Tz)'^ is a vector of Pauli matrices acting on the spin index. Non-conservation 
of the spin current in F^ and F^ results in a torque Ta {rt) exerted on the magnetic moment of 
-Pa (Fb) (the origin is chosen in the normal metal spacer separating Fa and F^), 
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(5) 



3 Non-equilibrium torque 

For a wide trilayer with a diffusive spacer layer, the main contribution to the non-equilibrium 
torgue lies in the plane spanned by the two magnetic moments {ha,hb) (the xz plane in Fig. 
|l]).Q In the absence of the superconducting contact, the non-equilibrium torque t^/J^n induced 
by a current / changes of sign at 9 = 0, vr, thus stabilizing the configurations = or = n 
depending of the direction of the current. El The presence of the superconductor, however, gives a 
geometric constraint that allows one to find a relation between r°^ and r^^ , from which it follows 
that T^^ also vanishes at ^ = tt/2. To find this relation, note that the spin current vanishes 
inside S (for bias voltage eV < Aq). Although the ferromagnetic layers do not conserve spin 
current, they do conserve the majority and minority currents parallel to their exchange fields 
ha{b) , hence 

fa -L /la , fb _L hb- (6) 



Since J+oo = and /15 points along the z-axis, we find from Eq. (^, 

Jo = Jo X (7) 

Finally, we use that, for wide trilayers, J_oo is polarized along ha and therefore, by Eq. (^) gives 
no contribution to the torque, f] The torque Tq is thus given by the projection of Jq to the unit 
vector perpendicular to ha and hence, 

rr = rr cos ^. (8) 

Equation (|8|) shows that the non-equilibrium torque r^*' can stabilize the 6 = tt/2 configuration. 
(Note, however, that the presence of the S contact creates an asymmetry between the torques 
Ta and Th, so that it is necessary for this effect that the magnetic moment of F;, is held fixed.) 

4 Equilibrium torque 

The previous section dealt with the non-equilibrium torque, which lies in the plane spanned by 
the two magnetic moments ha and hf,. It was found that r^*^ = for 6 = tt/2. We now ask 
whether the equilibrium torque has the same property. 

The answer is negative. Again, we can see this by a geometrical argument. In equilibrium, no 
spin current flows outside the trilayer on either side, and thus f^^ = —7^^ = Jo- In combination 
with Eq. (^), one then finds that r^"^ and f^^ point out of the plane spanned by ha and hb. There 
is no relation as simple as Eq. (^) for the out-of-plane component of the torque, hence no special 
behavior at ^ = 7r/2 is expected. 

There are two ways to compute the equilibrium torque: either by a direct calculation of the 
equilibrium spin currents in the trilayer, cf. Eq. (^), or from the derivative of the ground state 
energy with respect to the angle 9, 

It should be emphasized that the equilibrium torque involves contributions from energies in 
the entire conduction band, as opposed to the non-equilibrium torque where only quasiparticle 
states near the Fermi level play a role. 

To illustrate the dependence of t^, we have performed a numerical exact diagonalization 
of the BdG equation (|I|). We have considered the special case of a one-dimensional trilayer, 
i = ix va. Eq. (||). In Fig. g, Ta^ is plotted as a function of 9. We notice that the equilibrium 
torque (exchange interaction) depends strongly on the value of the superconducting gap Aq. 
This implies that both the direct spin exchange between the two magnets and other contribu- 
tions that involve the superconductor play a role. For a very large gap Aq ^ t, where 2t is the 
width of the conduction band in N, the superconductor acts as a hard wall and the results are 
qualitatively the same as without superconductivity. However, an interesting regime appear for 
relatively small values of the gap Aq < t where the equilibrium torque depends strongly on Aq, 
although the actual ^-dependence of t^ is very sensitive to the particular choice of parameters 
and, for a realistic calculation, would require knowledge of the detailed band structure of the 
different materials. In the particular example shown in Fig. y, for Aq = O.li the biquadratic 
harmonic dominates the exchange coupling 0, showing that, after eventual fine tuning of pa- 
rameters, nontrivial values of 6 can be stabilized. (A similar effect may also occur without a 
superconducting contact, see, e.g., Ref. 0.) 



"This is the case if the majority and minority spins are incoherently transmitted or reflected by the ferromag- 
netic layers, which is a good assumption for a wide trilayer with many propagating channels at the Fermi level. 
However, the conclusion (fel) also holds for coherent transmission or reflection Ref. pi 
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Figure 2: Equilibrium torque t""^ for a one- dimensional FNF trilayer with one superconducting contact. Torque 
is shown as a function of the angle 6 between the moments of the two ferromagnetic layers, and for various values 
the superconducting gap Ao, or without superconducting contact, see figure. The system consists of 80 sites, the 
relative sizes of the different layers being respectively for N, Fa, N, Ft, N and S: 9, 11, 9, 11, 4 and 36. We set 
Vi = €f ~ and hi — 0.8. All energies are measured in units of t. 
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